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Context. We wish to determine accurate ages for open clusters and use this, in conjunction with colour-magnitude 
diagrams, to constrain models of stellar structure and evolution. 
Aims. The detached eclipsing binary V20 in the old, metal-rich ([Fe/H] = +0.40) open cluster NGC6791 is studied in 
order to determine highly accurate masses and radii of its components. This allows the cluster age to be established 
with high precision, using isochrones in the mass-radius diagram. 

Methods. We employ high-resolution UVES spectroscopy of V20 to determine the spectroscopic orbit and time-series 
V, I photometry to obtain the photometric elements. 
Ck Results. The masses and radii of the V20 components are found to be 1.074 + 0.008A/q and 1.399 ± 0.016-Rq (primary) 

and 0.827 ± 0.004Mq and 0.768 ± 0.0067?q (secondary). The primary is located almost exactly at the hottest point 
along the cluster isochrone, and the secondary is a ~ 7 times fainter main-sequence star. We determine an apparent 
cluster distance-modulus of (m — M)v = 13.46 ± 0.10 (average of primary and secondary). The cluster age is obtained 
from comparisons with theoretical isochrones in the mass-radius diagram. Using the isochrones from Victoria-Regina 
with [Fe/H] = +0.37 we find 7.7 ± 0.5Gyr, whereas the Yonsei-Yale (V 2 ) isochrones lead to 8.2 ± 0.5Gyr, and BaSTI 
isochrones to 9.0 ± 0.5Gyr. In a mass-radius diagram, the 7.7Gyr VRSS and 9.0Gyr BaSTI isochrones overlap nearly 
perfectly despite the age-difference. This model dependence, which is significantly larger than the precision determined 
from mass, radius, and abundance uncertainties, prevents a definitive age-determination of the cluster. 
' Conclusions. Using detached eclipsing binaries for determination of cluster ages, the dominant error is due to differences 

among stellar models and no longer to observational errors in cluster reddening and distance. By observing a suitable 
, number of detached eclipsing binaries in several open clusters it should be possible to calibrate the age-scale and provide 

. firm constraints which stellar models must reproduce. 
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1. Introduction For these reasons, NGC 6791 has been the subject of a 

number of studies since the work by Kinman ([1965), and 
Among the old open clusters, NGC 6791 holds a special yet there is still not agre ement over its basic parameters. 
H I place: it is one of the oldest, most massive, and most Given a distance of ~ 4kpc and a non-negligible redden- 
metal-rich (Origlia [20061 Carretta, Braga glia fc Gratton in& determination of the cluster age is a complex prob- 
2007, Anthony-Twarog, Twarog & Meyer [20071 hereafter lem when attempting to use the colour-magnitude diagram 
ATTM07) clusters known. In addition to these features, (hereafter CMP). Even with r ecent attempts to detec t exo- 
the cluster cont ains a population of ho t blue stars (Liebert, planet transits (Brun tt et al. [20031 Mochejska et al. l2T)05l 
Saffer & Green [19941 Landsman et al. 11998] ) , a large pop- Montalto et al. M - which has led to a substantial body 
ulation of white dwarfs extending to the end of the cool- of well calibrated photometr y and a high-precision CMD 
ing sequenc e has been found with HST (Bedin et al. [20051 (Stetson, Bruntt & Grundahl [20031 hereafter SBG03) - it 
Bedin et al.[2008J), and extensive photometric studies have has proven very difficult to determine a precise cluster age. 
revealed a large p opulat ion of variable st ars in the cluster This is ma inly related to the difficulty of obtaining pre- 
fix grun tt et al. |20031 Mochejska et al. [2002J De Marchi cise reddening and distance estimates (see ATTM07 for an 
et al. | 2007| ) . extensive discussion) and to the problem of transforming 

model isochrones to observed colours and magnitudes at 

Send offprint requests to: F. Grundahl, the high metallicity of the cluster (Tripicco et al. H995|) . 
e-mail: fgj@phys.au.dk 

* Based on observations carried out at Nordic Optical The range of ages proposed for NGC 6791 has extremes 

Telescope at La Palma and ESO's VLT/UVES ESO, Paranal, of 7 and 12Gvr . but has m recent years seemed to converge 

Chile (75.D-0206A, 77.D-0827A). on a value near 8Gyr. This would put the cluster turnoff 

** Tables 11 and 12 are only available in electronic form at the mass close to IMq and thus bridge the "gap" between the 

CD S via anonymous ftp to cdsarc.u -strasbg.fr (130.79.128.5) or turnoff mass of globular clusters (typically 0.8Mq) and 

via |http: / / cdswebaT- strasbg.fr /cgi-bin/qcat7J/ A+ A/ 1 younger open clusters, where convective overshoot has a 
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marked effect on their turnoff morphology in the CMD. 
Such a high cluster age would also seem to suggest that 
atomic diffusion might have had sufficient time to act in 
its stars, such as found in the globular cluster NGC 6397 
(Korn et al. 120051) . 

In the course of a study to try and detect exoplanetary 
transits in the cluster (Bruntt et al. 2003) , we realized that 
the detached eclipsing system V20 (Rucinski, Kaluzny & 
Hilditch 1996) might offer the possibility to determine the 
masses and radii for stars near the turnoff. Rough estimates 
of the system parameters (Bruntt et al. I2003[) suggested 
that although V20 is only at V = 17.34, it would still be 
within the capabilities of UVES (Dekker et al. |2000|) at 
the ESO VLT to measure precise radial velocites for its 
components. 

It is well known that detached eclipsing binaries offer 
the possibility to determine precise masses and radii for 
the system components. If one or both components is close 
to its turnoff mass (for the age of the binary) , it is possible 
to put tight constraints on the age of the system through a 
comparison of the position of the primary and secondary in 
a mass-radius (M — R) diagram to theoretical isochrones. 
For stellar clusters, such an analysis has some significant 
advantages: the determination of the masses and radii is 
essentially independent of the usual "trouble-makers" such 
as reddening, distance, metallicity. Furthermore, since the 
comparison to models is carried out in the M — R diagram, 
one avoids the difficult process of transforming the effective 
temperatures and luminosities of the models to observed 
colours and magnitudes. Thus, a determination of cluster 
ages in the M — R diagram is essentially "the closest" one 
can get to a "direct" confrontation between observations 
and models. 

In this paper we undertake a full analysis of the de- 
tached eclipsing binary system V20 in NGC 6791 and deter- 
mine accurate values for the masses and radii of its compo- 
nents. This is used to determine a precise cluster age and to 
compare the ages derived using commonly used isochrones 
from three sources. Our main conclusion is that using such 
binary systems as V20, it is the "accuracy" of the available 
theoretical stellar models, which limit the obtainable age 
precision, and not the observational data. 



2. Photometry of V20 

The photometric data for V20 consists of V (Johnson) and / 
(Cousins) CCD observations from the 2.56m Nordic Optical 
Telescope (NOT) and its stand-by camera StanCam. We 
refer to the telescope homepage for further information^. 
V20 was observed on 9 nights between April 2003 and July 
2004, and a total of 298 and 300 exposures were obtained 
in V and /, respectively. Primary eclipse observations were 
done on three nights, secondary eclipse observations on four 
nights, and out-of-eclipse phases on two nights. Although 
complete coverage of the light curve has not been secured 
outside eclipses, the data are sufficient for a detailed pho- 
tometric analysis. 

For all observations we employed an exposure time of 
240s in V and 180s in /. Seeing conditions varied substan- 
tially from night to night, with median values of l'.'OO and 
0'.'82 in V and /, respectively. The best seeing frames have 
a FWHM of 0"60 and 0'.'46, respectively. 
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Fig. 1. The field of V20 in NGC 6791, centred at 
(a,<S)20oo = (19 71 20 m 54.30 s + 37 d 45 m 34.7 s ). Note, that 
in the left image V20 appears slightly elongated in the N-S 
direction, the FWHM in this I band image is 0'.'46. The 
righthand image shows the star subtracted version of the 
left image. Note the larger than expected residuals at the 
position of V20, consistent with a fainter third light compo- 
nent. The pixel size is 0'.'1755 and the greyscale is identical 
for the two images. 



The bias frames and flat fields, used in the data re- 
duction, were obtained during evening twilight on each 
observing night. All photometry was carried out with 
D AOPHOT /ALLSTAR/ ALLFR AME fStetson fTWlfTOO^ 
and DAOGROW (Stetson |1990p and transformed to a com- 
mon coordinate system using MATCH and MASTER (P. 
Stetson, private comm.). For each frame we produced a 
point-spread function (PSF) using the brightest stars in 
the field, and subsequently carried out aperture photometry 
in large apertures (neighbour stars subtracted) using the 
NEDA routine provided with DAOPHOT. Subsequently 
the aperture photometry was fed into DAOGROW to ob- 
tain the final large- aperture magnitudes. We found that this 
prodedure gives slightly better photometric precision than 
using profile-fitting photometry. 

During the reductions of the data, we noticed on the 
best seeing frames (FWHM — 0'.'46) that V20 appeared 
slightly elongated, and we interpret this as a signature of 
a third light component due to a chance alignment with a 
cluster star. This has to be properly taken into account in 
the final analysis of the lightcurves. In FigQ] we show the 
field in NGC 6791 centred on V20. Attempts were made to 
carry out ALLFRAME photometry with two components 
at the position of V20 separated by 074. This improved 
the residuals in the fit for the best seeing frames, but the 
large seeing variations did not allow a sufficient precision 
in the magnitudes of the V20 primary+secondary and the 
third light component. We therefore include the light from 
the third star in the light curves, and deal with it in the 
subsequent analyses. 

2.1. Light curves and standard indices 

With the instrumental photometry in hand, we proceeded 
to transform the observations to the V and / standard sys- 
tem. In SBG03 a large effort was put into this transforma- 
tion, and we have used the available photometry from this 
source as internal standard stars. For each frame, a linear 
transformation from instrumental magnitudes to standard 
magnitudes, using (V — I) as colour term, was calculated, 
and subsequently we averaged the coefficient for the colour 
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term for all frames (for a given filter) and used this for the 
final determination of the zeropoint for each frame. Our 
photometry is therefore on the same system as SBG03 with 
the out-of-eclipse magnitudes given by their photometry as 
listed in Table [TJ The accuracy of the photometry is in the 
range (P01 to ff n 02 as mentioned in SBG03. 

The light curves for V20 are listed in Tables 11 and 
12 and are shown in Figs. [2] and [3] with phases calculated 
from the ephemeris given in Eq. Q] V20 is clearly well de- 
tached with practically constant light level outside eclipses. 
Secondary eclipse occurs at phase 0.50, and the eclipses are 
of the same duration, supporting that the orbit is circu- 
lar. As seen from the different eclipse depths, V20 consists 
of two stars of rather different surface fluxes. Secondary 
eclipse is total; depths are (Fill (V") and m 142 (7). The 
close companion is included in the light curves, meaning 
that a significant amount of third light is present. 

Throughout the paper, the component eclipsed at the 
deeper eclipse at phase 0.0 is referred to as the primary (p), 
and the other as the secondary (s) component. 

Standard, out-of-eclipse, V, I photometry for V20 is 
listed in Table Q] together with calculated individual pho- 
tometry for the three stars. In order to obtain constraints 
on the contribution from the third light, we shall attempt 
to determine its value from the available light curve. 

We first determined the depth of the total secondary 
eclipse in V and I to be 0^111 and 0^142 , respectively. 
Since this represents the case where the secondary com- 
ponent of the binary is totally eclipsed, we can determine 
the magnitudes for the secondary component, and for the 
sum of the primary and third light component. If we as- 
sume (as is verified later from spectroscopy) that V20 is 
a member of the cluster, and that the third light is also 
a cluster member, then by assuming that it is located on 
the main sequence, a unique combination can be calculated, 
for which the sum of primary and third light matches the 
light during total eclipse. In practical terms, we fitted a low 
order polynomial to the main-sequence of the cluster (rep- 
resenting possible locations for the third light component). 
Then, for each point along the polynomial, the correspond- 
ing location of the primary component was calculated, such 
that the sum of their light matched the location of the light 
during total eclipse of the secondary component. 

The results of this exercise is given in TableQ] for the 
primary, secondary and third star. In Figd] we show the 
cluster (V — I), V colour-magnitude diagram based on the 
Stetson et al. (2003) photometry with the location of the 
V20 combined light and the individual components indi- 
cated. We have attempted to estimate the uncertainty in 
the location of the secondary component by assuming that 
the uncertainty in the estimate of the V and I depths of the 
total eclipse is 0™002 and then calculating its position 10000 
times at random from a normal distribution of a = 0™002 
in each colour and calculating the position of the secondary. 
As can be seen from Fig[J]it appears that three cluster stars 
can indeed match the observed total light (as expected) and 
that the third light component has a luminosity in between 
the primary and secondary. 

The exact luminosity for the third star is somewhat de- 
pendent on how the fiducial for the main-sequence is de- 
rived. We estimate that the uncertainty in its determined 
position is not larger than 0™15 - this is based on exper- 
iments where small shifts in colour were added to the fit- 
ted main-sequence polynomial, while it still appeared to 
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Fig. 2. V light curve for V20 (third star included). Only 
the observed phase intervals are shown. 
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Fig. 3. I light curve for V20 (third star included). Only 
the observed phase intervals are shown. 

Table 1. Standard V,I photometry for V20 outside 
eclipses from SBG03; V20 is star number 8600 in their list. 
The mean values, with errors, are based on 925 (V) and 703 
(I) observations, respectively. Photometry for the two com- 
ponents and the 3rd star was calculated from the depths of 
the total secondary eclipse, assuming that all three stars 
are located on the main sequence of NGC6791; see text for 
details. 



Object 


V 


I 


V-/ 


V20 + companion (SBG03) 


17.3390 


16.3560 


0.9830 


±5 


±6 


±8 


Primary 


17.6400 


16.7277 


0.9123 


Secondary 


19.8686 


18.6353 


1.2333 


3rd star 


19.4322 


18.2911 


1.1322 



represent the main-sequence. We note, that high-resolution 
imaging would allow a direct determination of the contribu- 
tion to the total light from the third star. For the expected 
separation (~0'/4), this would probably require adaptive 
optics. 
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Fig. 4. The (V — I, V) colour-magnitude diagram for 
NGC6791 with V20 and its components overplotted. The 
filled square indicates the location of the total light of the 
entire system. The light of the primary and third compo- 
nent (light during total eclipse of the secondary) is indi- 
cated by an open square. In order of decreasing luminosity, 
the filled circles indicate the location of the primary, the 
third star, and the secondary. For the secondary, we have 
overplotted the possible locations if errors of 0™002 in V 
and /-depths of the total eclipse are assumed. 

Table 2. Times of minima for V20. O-C values are calcu- 
lated for the ephemeris given in Eq. [T] adopting a circular 
orbit. Refere nces a re: M2002, Mochejska et al. l2"UU2l B2003, 
Bruntt et al. l2"0"0"51 G2008, this paper. 



HJD 

-2400000 


rms 


Type 


O-C 


Band 


Ref. 


52109.77415 


0.00110 


P 


0.00215 


R 


M2002 


53151.60485 


0.00030 


P 


-0.00161 


V 


G2008 


53151.60449 


0.00100 


P 


-0.00125 


I 




53180.54732 


0.00060 


P 


0.00153 


V 




53180.54747 


0.00070 


P 


0.00138 


I 




52102.53467 


0.00100 


S 


-0.00373 


V 


B2003 


52102.53332 


0.00080 


s 


-0.00238 


V 




52811.56249 


0.00100 


s 


0.00101 


V 


G2008 


52811.56404 


0.00090 


s 


-0.00054 


I 




52869.44066 


0.00210 


s 


0.00039 


V 




52869.44309 


0.00170 


s 


-0.00204 


I 





2.2. Times of minima and ephemeris 

From the V and I light curve observations, two times of 
minima per band have been derived for both eclipses. They 
are listed in Table together with times redetermined from 
the photometry by Bruntt et al. (120030 and Mochejska et 
al. (|2002[) . The method of Kwee and van Woerden (|1956|) 
was applied throughout. 



Weighted linear least square fits to the times of min- 
ima yield formally identical periods of 14 d 46989 ± 0.00005 
and 14.46999 ±0.00017 for primary and secondary eclipses, 
respectively. Since only three times per eclipse are avail- 
able, we have independently determined both the epoch and 
the period from JKTEBOP analyses of our V and / light 
curves; see Sect. 12.31 The periods obtained from the two 
light curves agree well and are close to those given above. 
We adopt the following linear ephemeris for all analyses in 
this paper: 



Min I = 2453151.6061 
±9 

2.3. Photometric elements 



14<?469918 x E 
±25 



(1) 



Since V20 is well detached with relative component radii of 
only about 0.045 and 0.025, respectively, we have adopted 
the simple Nelson-Davis-Etzel model (Nelson & Davis 1972, 
Etze HT9"5TI Martynov[1973) for the light curve analyses. It 
represents the deformed stars as biaxial ellipsoids and ap- 
plies a simple bolometric reflection model. We have used the 
corresponding JKTEBOF0 code, which is a revised and ex- 
tended version of the original EBOP code (Etzel fTMTj) . The 
Levenberg-Marquardt minimization algorithm (MRQMIN: 
Press et al. U992p is used for the lest-squares optimization of 
the parameters, and the code has been extended to include 
also non-linear limb darkening and adjustment of epoch and 
orbital period. In one of its modes, JKTEBOP performs 
Monte Carlo simulations, which we use to assign realistic 
errors to the photometric elements. For further information, 
we refer to e. g Sout hworth et al. (|2004a[ [2004bl [20071) and 
Bruntt et al. l|2"005|l . 

The V and I light curves were analysed independently, 
with equal weights assigned to all observations. The mag- 
nitude at quadrature was always included as an adjustable 
parameter, and the phase of primary eclipse was allowed to 
shift from 0.0. In initial JKTEBOP analyses, the epoch 
and orbital period was included as free parameters and 
then fixed at the values given in Eq. [TJ see Sect. 12.21 A 
circular orbit was assumed throughout, and the mass ra- 
tio between the components was kept at the spectroscopic 
value (Table [5]) . Gravity darkening coefficients correspond- 
ing to convective atmospheres were applied, and the simple 
bolometric reflection model built into EBOP/ JKTEBOP 
was used. Linear limb darkening coefficients were either 
assigned from theoretical calculations (Van Hamme[1993; 
Claret I2000|) according to the effective temperatures, sur- 
face gravities, and abundance, or left free. 

In text and tables on photometric solutions we use the 
following symbols: i orbital inclination; r relative radius; 
k = r s /r p ; u linear limb darkening coefficient; y gravity 
darkening coefficient; J central surface brightness ratio; L 
luminosity; I3 third light fraction. 

As mentioned in Sect. I2"land l2.1l a close star is included 
in the light curve data. The corresponding amount of third 
light I3 = Ls/(L P + L s + L3), as calculated from the infor- 
mation in Table [TJ is 0.146 ± 0.021 (V) and 0.168 ± 0.024 
(I) . Uncertainties are based on realistic estimates of the ac- 
curacy of the magnitudes of the total light (±0'. n 010), the 
depths of secondary eclipse (±0™002), and the magnitudes 
calculated for the companion (±07*15). We have adopted 
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Table 3. Photometric solutions for V20 from the 
JKTEBOP code. A photometric scale factor (the magni- 
tude at quadrature) and the phase of primary eclipse were 
included as free parameters. Linear limb darkening coef- 
ficients from Claret ([5000), van Hamme (|1993p . and free, 
respectively. The errors quoted for the adjusted parameters 
are the formal errors determined from the iterative least 
squares solution procedure. 



Band 


V 


V 


V 








Limb 


COO 


VH93 


Free 


COO 


VH93 


Free 


darkening 














i n 


89.76 
±11 


89.74 
±9 


89.85 
±17 


89.78 
±16 


89.78 
±15 


89.86 
±28 


r p + r B 


0.0703 
±3 


0.0695 
±3 


0.0696 
±4 


0.0709 
±5 


0.0703 
±5 


0.0705 
±5 


k — r s 1 r v 


0.539 
±1 


0.548 
±1 


0.545 
±3 


0.544 
±2 


0.553 
±2 


0.547 
±4 


r p 


0.0457 


0.0449 


0.0450 


0.0459 


0.0454 


0.0456 


r 3 


0.0246 


0.0246 


0.0245 


0.0250 


0.0250 


0.0249 


Up 


0.72 


0.61 


0.65 
±3 


0.56 


0.44 


0.51 
±4 


u s 


0.79 


0.74 


0.34 
±25 


0.62 


0.54 


0.43 
±26 


Vp 


0.39 


0.39 


0.39 


0.25 


0.25 


0.25 


Vp 


0.46 


0.46 


0.46 


0.29 


0.29 


0.29 


Js/Jp 


0.448 
±3 


0.447 
±3 


0.381 
±35 


0.599 
±5 


0.594 
±4 


0.562 
±56 


L s /L p 


0.126 


0.127 


0.127 


0.173 


0.174 


0.174 


h 


0.146 


0.146 


0.146 


0.168 


0.168 


0.168 


'(mag.) 


0.0044 


0.0044 


0.0044 


0.0067 


0.0067 


0.0067 



these 13 results, which are based on the realistic assump- 
tion that the third star is a main sequence cluster member, 
since it turns out that the amount of third light is not well 
constrained by the light curves; see below. 

The photometric solutions for different adopted theoret- 
ical linear limb darkening coefficients are given in Table [3l 
and (O — C) residuals are shown in Figures [5] and O which 
clearly reveal the higher quality of the V light curve. Except 
for a larger scatter, partly due to small nightly differences, 
there are no serious systematic trends in the / residuals. 
As seen, the orbital inclination is close to 90°, and the 
elements of the individual solutions agree well. The two 
quite different components are small compared to the ra- 
dius of the orbit; for both components the individual radii 
obtained agree within about ±1%. Additional tests reveal 
that the ratio of relative radii k is well constrained by the 
light curves. 

Including the linear limb darkening coefficients for both 
components as adjustable parameters was only partly suc- 
cessful, since those for the secondary components converged 
towards unrealistically low and also very uncertain values. 
For the primary component, the V value is close to that by 
Van Hamme (|1993[) . whereas the / value is just between the 
two theoretical coefficients. They are almost independent of 
whether u s is fixed at the theoretical values or adjusted. 

Adopting a non-linear limb darkening law (square root) 
was also attempted but has no significant influence on the 
solutions, as also found by Lacy, Torres & Claret (2008). 

As already mentioned the amount of third light is not 
well constrained by the light curves themselves. If included 
as free parameters, much too low but also very uncertain 
values are obtained (V : 0.04 ± 0.08; / : 0.13 ± 0.10). Upper 
third light limits, leading to an orbital inclination of 90°, are 
about 0.16 (V) and 0.18 (J). The effect on the photometric 
elements of changing the adopted Z3 values by 15%, i.e. by 
about their uncertainties given above, is shown in Table 2] 



Table 4. The effect of changing the amount of third light I3 
by 15%. Photometric solutions for V20 from the JKTEBOP 
code. A photometric scale factor (the magnitude at quadra- 
ture) and the phase of primary eclipse were included as 
free parameters. Linear limb darkening coefficients from 
van Hamme (|1993[) were adopted, see Table [3] The errors 
quoted for the adjusted parameters are the formal errors 
determined from the iterative least squares solution proce- 
dure. 



Band 


V 


I 


i n 


89.63 


89.62 


±6 


±9 


r p + r s 


0.0697 


0.0704 




±3 


±5 


h Ts/fp 


0.540 


0.542 


±1 


±2 


1-p 


0.0453 


0.0457 




0.0244 


0.0248 


J s 1 Jp 


0.446 


0.590 


±3 


±4 


L s 1 Lp 


0.123 


0.168 


h 


0.124 


0.143 


a (mag.) 


0.0044 


0.0067 



Table 5. Adopted photometric elements for V20. The in- 
dividual flux and luminosity ratios are based on the mean 
stellar and orbital parameters. 



i (°) 89.76 ±0.15 
r p 0.0452 ± 0.0005 
r s 0.0248 ± 0.0002 

V I 

Js/J P 0.451 0.590 

±7 ±8 

LJL„ 0.129 0.172 

±4 ±6 



The final, adopted photometric elements for V20 are 
listed in TableEl Higher weight has been given to the results 
based on Van Hamme (11993P linear limb darkening coeffi- 
cients. Uncertainties are based on interagreement between 
the V and / solutions, the uncertainties of the amount of 
third light, and Monte Carlo simulations. 

3. Spectroscopy for V20 

The spectroscopic observations for V20 were carried out in 
service mode with UVES at the ESO VLT during alloca- 
tion periods 75 and 77. Since V20 is at declination +37°, 
it can on Paranal only be observed at an airmass larger 
than 2.1 - therefore all observations were carried out near 
meridian passage. Due to the faintness of V20, and in order 
to minimize slit losses, a slit of 1'.'20 width, correspond- 
ing to a resolution of approximately 37 000, was used. The 
slit was aligned along the parallactic angle (ELEV mode), 
and the ADC was not inserted in the beam, since it causes 
a slight loss of flux. The standard 580nm setup, and on- 
chip binning of 2x2 pixels, was used for all observations. 
The corresponding wavelength ranges covered at the two 
CCD detectors employed in UVES are approximately 4775- 
5750 A and 5875-6830 A, respectively. A total of 17 epochs 
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Fig. 5. (0-C) residuals of the V observations from the 
theoretical light curves for the VH93 parameters listed in 
Table El 

were obtained, see Table [HI each with a ThAr exposure at- 
tached. The FWHM from the image headers for the start 
and end of each exposure are included. The typical S/N 
per pixel for the red chip was between 15 and 25; two of 
the spectra had to be omitted in the data analysis due to a 
very low signal. After careful check of the wavelength cali- 
brations, we decided to apply the pipeline reduced spectra 
for the analyses described below. 

3.1. Radial velocities and spectroscopic elements 

In order to reduce/eliminate possible systematic velocity 
errors, we have decided to apply the broadening function 
(BF) formalism (Rucinski QUEE Hffl EQQH with useful IDL 
routines and a cookbook from his homepag^fl) for radial 
velocity measurements of V20. Based on results reported 
in the literature, we expected this method to be well suited 
for the analysis of the sharp lined but otherwise non-trivial 
spectra of V20 with three set of lines. As seen below, our 
expectations have been fully met. 

The BF method applies only one (unbroadened) tem- 
plate at a time, so we have decided to perform several 
analyses adopting three different synthetic templates, cor- 
responding to each of the three set of lines present in the 



Fig. 6. (O-C) residuals of the / observations from the 
theoretical light curves for the VH93 parameters listed in 
Table H 

Table 6. Log of spectroscopic observations for V20 ob- 
tained with UVES, ordered by orbital phase. The helio- 
centric Date (HJD) is given at mid-exposure, exposure 
time (T exp ) in seconds, and seeing (FWHM) in arcseconds. 
Observations excluded in the analyses are marked by *. 



HJD 


Phase 


Texp 


Airmass 


FWHM 


2453847.88977 


0.11939 


4060 


2.34-2.15 


0.73-0.66 


2453573.65870 


0.16759 


2700 


2.20-2.15 


0.79-0.69 


2453588.62071 


0.20160 


2700 


2.19-2.15 


0.89-0.68 


2453588.65279 


0.20381 


2700 


2.15-2.25 


0.67-0.60 


2453589.60199 


0.26941 


2700 


2.27-2.16 


0.82-0.83 


2453589.63401 


0.27162 


2700 


2.16-2.18 


0.96-0.75 


2453590.62070 


0.33981 


2700 


2.18-2.16 


0.50-0.73 


2453590.65270 


0.34202 


2700 


2.16-2.28 


0.64-0.63 


2453880.81608 


0.39489 


4060 


2.24-2.16 


0.72- 


2453623.51512 


0.61311 


2700 


2.23-2.15 


0.87-0.87 


2453623.54710 


0.61532 


2700 


2.15-2.21 


0.98-0.81 


2453986.50872* 


0.69919 


4060 


2.37-2.15 


-1.11 


2453957.58875* 


0.70057 


4060 


2.37-2.15 


1.39-1.95 


2453943.67176 


0.73878 


4060 


2.16-2.24 


0.84-0.87 


2453871.82770 


0.77372 


4060 


2.32-2.15 


0.45-0.36 


2453626.50659 


0.81985 


2700 


2.23-2.15 


0.78-0.72 


2454003.51981 


0.87481 


4060 


2.15-2.33 


0.93-0.81 



F. Grundahl et al.: V20 in the open cluster NGC6791 



7 



Table 7. Spectroscopic elements for V20 determined from the two spectral regions applying different templates. Both 
double-lined and single-lined solutions are included. Epoch and orbital period were fixed at the values listed in Eq. [U 
and a circular orbit was assumed. K p and K s are the velocity semi-amplitudes for primary and secondary component, 
respectively, 7 is the system velocity for V20, and RV3 is the mean value of the radial velocity for the 3rd star, a is the 
standard error of one observation. All quantities are given in units of kms -1 . Spectra regions (R) are: L = 5100-5600 A; 
U = 6100-6700A with H Q blocked out. Templates (T) used are P (Primary): T cff = 5600 K, log(g) = 4.2, [Fe/H] = +0.4; 
S (Secondary): T eff = 4800 K, log(g) = 4.6, [Fc/H] = +0.4; C (3rd star): T cS = 5000 K, \og(g) = 4.5, [Fe/H] = +0.4. The 
solutions marked by * are based on velocities corrected by the results from the analyses of synthetic spectra; see text for 
details. 



R 


T 








K„ 


(Tn 


u s 




T 
I 




RVz 


L 


1 


47 


04 ± 


20 


01.22 ± 0.20 


0.37 


0.89 


—46 


49 + 


12 


A A AO _l_ f\ £1 1 

— 44. 4o ± U.ol 


T * 

L 1 




47 


01 ± 


19 


/"* "1 O 1 1 /"ill"! 

61.31 ± 0.19 


0.37 


0.86 


—46 


49 + 


12 


a a ao 1 n r* 1 

—44.43 ± 0.61 


T 

L 


c 

a 


47 


15 ± 


20 


ol.lo ± 0.20 


0.39 


0.92 


—46 


46 + 


12 


A A CO J ACT 

— 44.5z ± U.57 


L 


C 


47 


13 ± 


20 


61.14 ± 0.20 


0.38 


0.89 


—46 


47 + 


12 


A A V (t 1 C\ f T 

—44.52 ± 0.57 


U 


P 


47 


19 ± 


11 


61.24 ± 0.11 


0.35 


0.43 


—46 


56 + 


07 


A A f r 1 n TO 

—44.55 ± 0.73 


u* 


p* 


47 


07 ±0 


11 


61.18 + 0.11 


0.37 


0.41 


-46 


50 + 


07 


-44.50 ± 0.73 


u 


S 


47 


18 ±0 


15 


61.16 + 0.15 


0.36 


0.62 


-46 


54 + 


09 


-44.49 ± 0.71 


u 


c 


47 


19 ±0 


13 


61.15 + 0.14 


0.36 


0.57 


-46 


54 + 


08 


-44.49 ± 0.70 


L 


p 


46 


98 ±0 


08 




0.27 




-46 


71+0 


07 




L* 




46 


96 ±0 


08 




0.27 




-46 


71+0 


07 




L 


s 


47 


08 ±0 


08 




0.26 




-46 


71 + 


07 




L 


c 


47 


06 ±0 


08 




0.26 




-46 


71 + 


07 




U 


p 


47 


17 ±0 


10 




0.33 




-46 


63 + 


09 




u* 




47 


04 ±0 


10 




0.33 




-46 


63 + 


09 




u 


s 


47 


17 ±0 


11 




0.35 




-46 


62 + 


09 




u 


c 


47 


17 ±0 


10 




0.34 




-46 


62 + 


09 




L 


p 








61.16 + 0.25 




0.82 


-46 


26 + 


22 




L 


s 








61.07 + 0.26 




0.84 


-46 


20 + 


22 




L 


c 








61.07 + 0.25 




0.81 


-46 


22 + 


21 




U 


p 








61.22 + 0.13 




0.41 


-46 


50 + 


11 




u 


s 








61.14 + 0.18 




0.56 


-46 


48 + 


16 




u 


c 








61.13 + 0.17 




0.54 


-46 


47 + 


14 





0.08 



0.06 



0.04 



0.02 



0.00 




-120 



-80 -60 -40 -20 
Velocity (km/s) 



Fig. 7. Broadening function obtained for the 6100-6700 A 
region using the primary (P) template. The UVES spec- 
trum was taken at phase 0.342, at HJD=2453590.65270. 
From left to right the components are: primary, third star, 
secondary. 



Table 8. Spectroscopic orbital solution for V20. T is the 
time of central primary eclipse. 



Parameter 



Value 



Adjusted quantities: 




K p (kms -1 ) 


47.09 + 0.09 


K s (kms" 1 ) 


61.16 + 0.20 


7 (kms -1 ) 


-46.63 + 0.13 


Adopted quantities: 




T (HJD-2 400 000) 


53151.6061 


P (days) 


14.469918 


e 


0.00 


Derived quantities: 




M p sin 3 i (Mq) 


1.074 ± 0.008 


M s sin 3 z (M ) 


0.827 ± 0.004 


asini (Rq) 


30.946 ± 0.063 


Other quantities: 




N ob3 


15 


Time span (days) 


377 



V20 spectra. For calculation of templates we have used the 
bssynth tool (Bruntt, private communication), which ap- 
plies the SYNTH software (Valen ti fc P iskunov [1996]) and 
modified ATLAS9 models (Heiter EUMJ) . The UVES spec- 



tra were first carefully cleaned and normalized, and then, 
like the template spectra, logarithmically rebinned to a con- 
stant velocity step of 1.0 kms -1 . 

The two spectral regions were analysed separately, and 
the position of the BF's, defining the radial velocities of the 
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Fig. 8. Spectroscopic double-lined orbital solution for V20 
obtained from the 6100-6700 A region using the primary (P) 
template; see Table [7] Filled circle: primary; open circle: 
secondary. The horizontal dotted line (upper panel) rep- 
resents the center-of-mass velocity of V20. The + signs 
are the measured velocities for the third light component. 
Phase 0.0 corresponds to central primary eclipse. 



stars, were calculated by fitting Gaussians to the smoothed 
functions. A sample BF is shown in Fig. [7] 

Spectroscopic elements were derived from the measured 
radial velocities using the method of Lehman-Filhes imple- 
mented in the SBOP program (Etzel 2004)) . The orbital pe- 
riod given in Eq. [T]was adopted, and the orbit was defined 
to be circular. Both double-lined and single lined solutions 
were done. Assuming that the radial velocity of the third 
component is constant, small shifts corresponding to the 
difference between its mean velocity and the individual val- 
ues were added to measured velocities of the primary and 
secondary components. The shifts do not alter the derived 
system velocities and semi-amplitudes significantly but im- 
prove their uncertainties. 

The orbital elements determined from the different spec- 
tral ranges using different templates are given in TableUJto- 
gether with the mean radial velocities for the third star. As 
seen, very accurate velocity semi- amplitudes are obtained; 
in general, the highest accuracy is reached for the 6100-6700 
A region; Fig.[8]illustrates one of the solutions. The double- 
lined solutions and the single-lined solutions from the 6100- 
6700 A region agree extremely well, whereas the semi- 
amplitudes of the single-lined solutions from the 5100-5600 
A region are systematically lower by about 0.10 kms -1 . 
The primary and secondary single-lined solutions, except 
the secondary for the 5100-5600 A region, give slightly dif- 
ferent values for the system velocity, although the results 
formally agree within the errors. 



In order to check for possible systematic velocity er- 
rors, we have analysed synthetic spectra calculated for the 
observed orbital phases. The three templates, broadened 
to vsxai values matching the observations, were shifted to 
their velocities at a given phase and combined according 
to the relative luminosities of the V20 components and the 
third star. BF analyses of the synthetic spectra show that 
the measured radial velocities deviate only slightly from the 
input values, typically below ±0.15 kms -1 for the primary 
and ±0.30 kms -1 for the secondary. Furthermore, adding 
the small shifts to the radial velocities measured from the 
observed spectra does not lead to significant changes of 
the orbital elements and their uncertainties. Two examples 
(marked by *) are included in Table [7) 

We adopt the weighted mean of the uncorrected single- 
lined solutions as the final spectroscopic elements for V20 
listed in Table [5] Errors typical for the individual solutions 
have conservatively been assigned to the semi-amplitudes 
and the system velocity. As seen, minimum masses accu- 
rate to 0.7% and 0.5% have been obtained for the primary 
and secondary components, respectively. For comparison, 
masses derived from the corrected single- lined solutions are 
only 0.2% higher for the primary and 0.2% lower for the 
secondary component. 

We have determined a systemic velocity for V20 of 
-46.63 ± 0.13 km s^ 1 (Table EJ). This is in excellent agree- 
ment with the value of —47.1 ± 0.8 kms _1 as determined 
from 15 cluster members by Carraro et al. (12006P who also 
found the dispersion in the radial velocities of the 15 stars 
to be 2.2 ± 0.4 km s _1 . This leaves little doubt that the 
system is indeed a member of NGC 6791. 



4. Absolute dimensions and distance for V20 

Absolute dimensions for the components of V20 are calcu- 
lated from the elements given in Tables [5] and [8] As seen in 
TablelHl masses and radii have been obtained to an accuracy 
of about 0.6% and 1.0%, respectively. Individual V,I,V — I 
magnitudes and indices are included, as calculated from the 
combined photometry (Table [1]) and the luminosity ratios 
between the components (Table[5]). An uncertainty of 15%, 
correlated in V and /, has been assumed for Z3. 

For the determination of effective temperatures from the 
(V — I) indices, we adopt an interstellar reddening of E(B — 
V) = 0.15±0.02, see Sect.EU and assume [Fe/H] = ±0.40± 
0.10, see Sect. 15.21 The calibration by Ramirez & Melendez 
PUU5]) then gives 5715 ± 125 K and 4750 ± 150 K for the 
primary and secondary components, respectively, whereas 
5665 ± 100 K and 4900 ± 100 K ar e obta ined from the recent 
calibration by Casagrande et al. (|2006p . From the empirical 
flux scale by Popper (|1980p . the bolometric correction (BC) 
scale by Flower (1996), and the V flux ratio between the 
components, which is obtained to high precision from the 
light curve analyses (Table [5]), we get a well constrained 
temperature difference of 765 ±15 K, supporting the results 
from the Casagrande et al. calibration. We consequently 
adopt 5665 and 4900 K as our final T c s results. 

As seen in Table |H] the synchronous rotational veloci- 
ties of the components of V20 are small. Synthetic binary 
spectra based on these velocities and a resolution of 37 000 
qualitatively agree well with the observed spectra, but we 
have not attempted to determine the actual rotational ve- 
locities. 
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Table 9. 

B.C.Q = 



Astrophysical data for V20. T e 
-0.08, and M, 



5780 K, 5.2. Cluster metallicity 



bolQ 



We have adopted E(B - V) 
1.3 x E(B - V), and A v = 3.1 x E(B - V) 
equatorial velocity for synchronous rotation. 



4.74 has been assumed. 
0.15 ± 0.02, E(V -I) = 
x E(B - V). v sync is the 



Primary 



Secondary 



Absolute dimensions: 






M/M P) 


1.074 ± 0.008 


0.827 ± 0.004 


1 


1.399 ±0.016 


0.768 ± 0.006 


log q (cks) 


4.178 ±0.010 


4.586 ± 0.008 


Vsync (kmS _1 ) 


4.9 ± 0.1 


2.7 ±0.1 


Photometric data: 






V 


17.642 ± 0.028 


19.867 ± 0.041 


I 


16.729 ± 0.033 


18.639 ± 0.046 


(V-I) 


0.914 ± 0.008 


1.228 ± 0.044 


V 


17.177 ±0.068 


19.400 ± 0.074 


(V-I)o 


0.719 ±0.027 


1.031 ±0.051 


T cS (K) 


5665 ± 100 


4900 ± 100 


M hoi 


4.10 ±0.08 


6.03 ± 0.09 


log L/Lq 


0.26 ± 0.03 


-0.52 ±0.04 


B.C. 


-0.10 


-0.35 


M v 


4.20 ± 0.08 


6.38 ±0.09 


V -M v 


13.44 ± 0.09 


13.48 ±0.10 


Vo - M v 


12.98 ± 0.11 


13.02 ±0.12 


Distance (pc) 


3940 ± 200 


4015 ± 225 



Eclipsing binaries arc important primary distance indi- 
cators fClausen l2004p . and V20 provides a more direct de- 
termination for NGC6791 than traditional main-sequence 
fitting. The two components yield nearly identical distances 
(4000 pc or (V - M v ) = 13 m 00), which have been estab- 
lished to a robust accuracy of 5%, taking into account all 
error sources. For comparison, recently published results 
from main-sequence fitting, which show a spread signifi- 
cantly higher than the quoted errors, are: 12™79 (SBG03; 
E(B- V) = 0.09 and [Fe/H] = 0.3), 13™07 ±fT04 (Carney 
et al. 1217051 E(B - V) = 0.14 and [Fe/H] = 0.4), and 
13 m l ± (Fl ( ATTM07 ; E(B - V) = 0.155 and [Fe/H] = 
0.45). 



5. NGC6791 

5.1. The reddening of NGC6791 



Anthony-Twarog, Twarog & Mayer (|2007j) gave a thor- 
ough discussion of the reddening towards NGC 6791 based 
on vbyCaHP photometry and a recalibration suitable for 
metal-rich stars, such as those found in this cluster. Over 
the years, since the first study of Kinman (|1965p . there 
has been significant disagreement over the reddening value. 
Here we adopt a value of E(B - V) = 0.15 ± 0.02 for the 
cluster, giving most weight to the ATTM07 value which 
also agrees well with the value derived from the maps of 
Schlegel et al. 1(1955)1 . 



As is the case for the reddening towards NGC 6791, the 
determination of its metallicity also has a long history. We 
shall not repeat it here but refer the reader to ATTM07 and 
Carretta et al. (|2007|) for a more comprehensive discussion. 
Based on the intense interest over the past few years, we 
shall adopt the values determined spectroscopically. While 
a full consensus has not been reached on this subject, most 
recent determinations seem to agree on a value close to 
[Fe/H] = ±0.40, which we adopt together with an uncer- 
tainty of ±0.10 dex. 

We note that e.g. Carretta et al. report significant un- 
derabundances for C,N, and O of about -0.3 dex. Like 
Carraro et al. (2005)) they find scaled-solar a-element abun- 
dances. 



6. Comparison with theoretical models 

In the following, we compare the accurate dimensions ob- 
tained for V20, and the NGC 6791 CMDs, with properties 
of the VRSS (scaled-solar mixture) Victoria-Regina evolu- 
tionary tracks and isochrones (VandenBerg et al., [2006 f*l, 
the Y 2 (Yonsei-Yale) grids (Demarque et al. I2004[R and 
the extensive material available from the BaSTI database 
(Pietrinferni et al., [20040 These recent models are all 
based on up-to-date input physics but differ somewhat with 
respect to e.g. core overshoot and diffusion (if included) 
treatment, He enrichment law, adopted solar mixture, and 
envelope convection calibration. Furthermore their asso- 
ciated colour transformations are different. We refer to 
Clausen et al. (|2008[) for further details, as well as to the 
original papers. 

Throughout we have adopted scaled-solar mixture for 
the heavy elements, and we have selected the available mod- 
els and isochrones which are closest to the adopted [Fe/H] = 
±0.40 ± 0.10 for the comparisons. 



6.1. V20 

For a given mass and [Fe/H], the observable properties of 
models at a given age like radius, effective temperature, 
and luminosity depend on the adopted input physics, in- 
cluding treatment of core and envelope convection, diffu- 
sion etc., and the assumed Y,Z relation. In order to con- 
strain such free model parameters, accurate results from 
many binaries are needed; cluster members are particularly 
valuable. From the binary perspective, masses and radii are 
the most direct parameters available, free of any scale de- 
pendent calibrations. So, the M — R diagram is well suited 
for isochrone tests, especially when the binary components 
are significantly different, as is the case for V20. In ad- 
dition, the T G ff — R plane allows tests of model tempera- 
tures, which for a given mass and [Fe/H] depend on e.g. 
abundance mixture leading to Z, on Y, and on surface con- 
vection efficiency. Adopting the usual assumption of coeval 
formation of the components from the same raw material, 
identical ages must be reached in these two planes, as well 
as in the M — log(L) plane. 



|http : //wwwl . cad c-ccda. hia-ih a.nrc-cnrc .gc . ca/| 
cvo/community/VictoriaReginaModels/ 
^ |http : //www. astro .ya le ■ edu/demarque/yystar .html 
^ |http : //www . te . astro . it/BASTI/index .php| 
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Table 10. Ages for the componets of V20 as deter- 
mined from M — R isochrones calculated for different mod- 
els and chemical compositions. For each model, the first 
line represent the available models closest to the observed 
[Fc/H]. All models include core overshoot; BaSTI standard 
and core overshoot models yield nearly identical results. 
Uncertainties due to mass and radius errors are about 0.3 
Gyr (primary) and 0.9 Gyr (secondary). 
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Fig. 9. V20 compared to Victoria-Regina VRSS mod- 
els for (X,Y,Z) = (0.63656,0.32344,0.04000), equivalent to 
[Fe/H] = +0.37 for [a/Fe] = 0.00. Isochrones from 5.0 to 
9.0 Gyr (step 1.0 Gyr) are shown. 



Fig. 10. V20 compared to Y 2 models for (X,Y,Z) = 
(0.6467,0.3122,0.0411), equivalent to [Fe/H] = +0.40 for 
[a/Fe] = 0.00. Isochrones from 5.0 to 9.0 Gyr (step 1.0 Gyr) 
are shown. 
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From the M — R comparisons shown in Figs. [51 [TQl and 
[TTlwe derive the component ages listed in Table [TU1 As seen, 
mass and radius uncertainties translate to age uncertain- 
ties of about 0.3 Gyr (primary) and 0.9 Gyr (secondary), 
respectively, for a given adopted composition. However, 
at [Fe/H]= 0.40, the ages resulting from the three grids 
range between 7.6 and 9.0 Gyr for the primary and be- 
tween 6.2 and 8.6 for the secondary. Within uncertainties, 
the Victoria-Regina and BaSTI isochrones predict identi- 
cal ages for the two components, although the age for the 
secondary is systematically lower than that of the primary. 
The Y 2 grids clearly predict different ages, again with lower 
values for the secondary. Age uncertainties due to an [Fc/H] 
uncertainty of 0.1 dex are comparable to those coming from 
masses and radii. 

We note that this trend of lower ages for the less mas- 
sive star is opposite to what has been found for several field 
G-type eclipsing binaries; see e.g. Popper (1997), Clausen 



Fig. 11. V20 compared to BaSTI models with overshooting 
for (X,Y,Z) = (0.657,0.303,0.0400), equivalent to [Fe/H] = 
+0.395 for [a/Fe] = 0.00. Isochrones from 5.0 to 9.0 Gyr 
(step 1.0 Gyr) are shown. 

et al. CEHMI): Torres et al. flSDDBD, and Clausen et al. (in 
prep.). In these cases, the seemingly higher ages may, al- 
though not yet firmly proved, be due to surface activity and 
correlated less efficient convection at the less massive com- 
ponents, which result in larger radius and lower effective 
temperature (but nearly unchanged luminosity). 

In Fig. rjJl we compare the Victoria-Regina, BaSTI, and 
Y 2 isochrones fitting the primary component. The first two 
are nearly identical, but for an age difference of as much as 
1.3 Gyr, whereas the shape of the Y 2 isochrone, at an age 
approximately between the two, is significantly different. 

Turning to the T c g — R plane shown in Figs. [13l QjJ and 
rp31 it is seen that the Victoria-Regina tracks fit both com- 
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Fig. 12. 

V20 compared to the VRSS (7.7 Gyr, dashed thick gray), 
Y 2 (8.2 Gyr, dashed-dot black), and BaSTI (9.0 Gyr, 
dashed black) isochrones which fit the primary component 
at the observed [Fe/H]; see Table [TU1 

ponents within errors. They are very slightly hotter than 
observed, and as seen in Fig. [T3] this is also the case for 
tracks calculated for a higher [Fe/H]— 0.49, which actually 
fall close to the [Fe/H]= 0.37 tracks. The Y 2 tracks are 
cooler than observed, and here tracks for a 0.1 dex lower 
[Fe/H] agree better. The BaSTI tracks, both standard and 
with core overshoot, fit the components of V20 perfectly 
well. 

In conclusion, both the the Victoria-Regina and the 
BaSTI models represent V20 well, but at ages which differ 
by about 1.3 Gyr. This is significantly higher than the preci- 
sion of about 0.5 Gyr, which can be reached from the avail- 
able mass, radius, and abundance information. The cause(s) 
for this difference is not clear but could perhaps be related 
to the fact, that the Victoria-Regina models apply different 
Y, Z relations and core overshoot treatment. Direct com- 
parisons between Victoria-Regina and BaSTI tracks and 
isochrones are shown in Figs. [TBI and 1171 

Additional NGC6791 binaries with component (s) be- 
tween 0.85 and 1.0 Mq would add further constraints on 
the models, and would e.g. clearly reveal if the isochrone 
shape is correct. Also, similar comparisons for younger 
clusters would be valuable. We note that for 1-2 Gyr 
field F-type binaries with component masses in the 1.1-1.4 
Mq range, Clausen et al. (2005]) recently found that the 
Victoria-Regina models are superior to the BaSTI models. 

6.2. NGC6791 colour-magnitude diagrams 

In light of the results obtained from V20 in the previous 
section, it is interesting to investigate what extra informa- 
tion can be gained from a comparison of the isochrones to 
the CMD of the cluster. There is an abundant literature on 
age determination of old open- and globular clusters using 
"isochrone-fitting" . The method is well described and so 
are the problems involved in the process. Observationally 
it is necessary to obtain accurate determinations of red- 



Fig. 13. V20 compared to Victoria-Regina VRSS mod- 
els for (X,Y,Z) = (0.63656,0.32344,0.04000), equivalent 
to [Fe/H]= +0.37 for [a/Fe] = 0.00. Tracks for 0.8, 
0.9, 1.073, and 1.1 Mq (full drawn) and isochrones 
from 5.0 to 9.0 Gyr (dashed, step 1.0 Gyr) are shown. 
For comparison, tracks for 0.8, 0.9, and 1.1 Mq 
for (X,Y,Z) = (0.66856,0.30144,0.03000), equivalent to 
[Fe/H] = +0.23 for [a/Fe] = 0.00 (dotted), and for (X,Y,Z) 
= (0.60456,0.34544,0.05000), equivalent to [Fe/H] = +0.49 
for [a/Fe] = 0.00 (dot-dash), are included. 
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Fig. 14. V20 compared to Y 2 models for (X,Y,Z) = 
(0.6467,0.3122,0.0411), equivalent to [Fe/H] = +0.40 for 
[a/Fe] = 0.0. Tracks for the component masses (1.074 
and 0.827 Mq, full drawn) and isochrones from 5.0 to 
9.0 Gyr (dashed, step 1.0 Gyr) are shown. To illus- 
trate the effect of the abundance uncertainty, tracks for 
[Fe/H] = +0.30 and [a/Fe] = 0.00 (dotted), corresponding 
to (X,Y,Z) = (0.6686,0.2976,0.0338), and [Fe/H] = +0.50 
and [a/Fe] = 0.00 (dash-dot), corresponding to (X,Y,Z) = 
(0.6200,0.3300,0.0500), are included. 
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Fig. 18. Isochrones for ages 7.7 (VRSS, red), 8.2 (Y 2 , blue), and 9.0Gyr (BaSTI, black) overplotted on the various 
CMDs of NGC 6791 (SBG03 photometry, plotted at light grey squares). In all plots we have adopted E(B - V) = 0.15 
and (m — M)y — 13.46 to transpose the isochrones to the observational plane (see Table IH])- In the middle panel, the 
location of the primary and secondary component of V20 are shown as filled black circles. The models are the same as 
those used in Fig. [TJ] On black and white printers the models are hard to distinguish, but at the cluster turnoff the 
VRSS models are always the bluest and the Y 2 models always the reddest. 



dening, distance and cluster abundance, and subsequently 
match theoretical models shifted in colour and luminos- 
ity to the observed CMD. The theoretical models must be 
transformed to the observational plane (colour and magni- 
tude) in order to allow the comparison. It is well known that 
for old clusters, the isochrone shapes only change slowly 
with time, and thus age-determinations based on isochrone 
shapes are not very sensitive. Furthermore, the turnoff lu- 
minosity and colour also changes slowly with time. In the 
best cases, the accuracies obtained for reddening and dis- 
tance modulus of clusters is of the order 0^01 and 0™1, 
respectively (Pasquini et al. I2008p . For an 8Gyr old clus- 
ter such errors translate into an uncertainty in age (from 
the turnoff position) of ~ 0.5 Gyr and ~ 1.0 Gyr, respec- 
tively (estimated from VRSS isochrones and for [Fe/H] = 
0.4). Thus, to obtain age estimates better than 1 Gyr, very 
precise values for reddening and distance must be obtained 
- in addition to this the models also show a non-negligible 
dependence on the adopted metallicity. 

VandenBerg & Stetson (|2004[) . Fig. 7, present a system- 
atic investigation of what effect plausible changes in metal- 
licity and distance modulus (at fixed reddening) has on the 
derived age for NGC 188, which is ~2Gyr younger than 
NGC 6791. Their exercise shows that ages between 5.9 and 



8.1 Gyr are compatible with the cluster CMD correspond- 
ing to ~ 30% of the cluster age. 

To confront the theoretical models with the observa- 
tions we show, in Fig. [TSJ the CMD for NGC 6791 in the 
(B - V, V), (V - I, V) and [B - I, V) planes with the pre- 
ferred age for each of the three isochrone sets (using the 
reddening and distance from Table[5]). We employ the mod- 
els closest to [Fe/H] = +0.40 throughout. It is evident that 
the isochrones are not far from matching the CMD in all 
colour planes, although the Y 2 models are somewhat too 
red compared to the VRSS and BaSTI models. From in- 
spection of the figure, it is clear that offsets within 0™02 in 
reddening and 0^10 in distance- modulus would bring the 
models into agreement with the cluster turnoff region. We 
thus conclude that for the three age values, the isochrones 
in the CMD are not in direct conflict with the observations, 
except perhaps for the (V — /) and (B — I) colours of the 
Y 2 models. By simple experimenting with VRSS models 
it is also clear that with (plausible) modifications to the 
adopted reddening and distance modulus, it is posible to 
match the turnoff region for the cluster for ages spanning 
at least the 6-9 Gyr range; the same conclusion would also 
hold for BaSTI or Y 2 models. 
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Fig. 15. V20 compared to BaSTI models with overshooting 
for (X,Y,Z) = (0.657,0.303,0.0400), equivalent to [Fe/H] = 
+0.395 for [a/Fe] = 0.00. Tracks for the component masses 
(1.074 and 0.827 Mq, full drawn) and isochrones from 5.0 
to 9.0 Gyr (dashed, step 1.0 Gyr) are shown. For compar- 
ison, tracks (dotted) for {X,Y,Z) = (0.682,0.288,0.0300), 
equivalent to [Fe/H] = +0.254 for [a/Fe] = 0.00, are in- 
cluded. 
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Fig. 16. Comparison between VRSS (thick blue) and 
BaSTI (thin red) isochrones at 5,7, and 9 Gyr for com- 
positions close to the observed [Fe/H]; see Table [TD] The 
components of V20 are included for comparison 

Given the current uncertainties in the derivation of 
cluster reddening, distance modulus, and metallicity, it is 
the opinion of the authors that we can only gain limited 
information on cluster absolute ages from their colour- 
magnitude diagrams. From the results presented here, V20 
provide much tighter constraints on the age of NGC6791 
than isochrone-fitting to the colour-magnitude diagram. In 
fact, it would be necessary to know the distance-modulus 
and intrinsic colour of the cluster turnoff to an accuracy 
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Fig. 17. Comparison between VRSS (thick blue) and 
BaSTI (thin red) 1.1 and 0.8 M tracks (full drawn) and 
5,7, and 9 Gyr isochrones (dashed), for compositions close 
to the observed [Fe/H]; see Table [TOl The components of 
V20 are included for comparison 

better than 0™05 and 0™01, respectively to match the 
0.5Gyr age precision. 

Unfortunately, we are faced with the situation that even 
in the M — R diagram, the models give very discrepant ages 
with differences of 15-20% (see Fig.[S]). Until this situation 
is clarified, it is the models which limit the precision with 
which stellar ages can be obtained. We speculate that in 
the case of NGC 6791 the discrepancy could be due to dif- 
ficulties in the description of convective energy transport. 

7. Summary and conclusions 

In this paper we have presented extensive photometric and 
spectroscopic observations of the detached eclipsing binary 
V20 in the old open cluster NGC 6791, and determined pre- 
cise masses and radii for its primary and secondary compo- 
nents. The errors in masses and radii are not subject to the 
usual problems of determining the cluster reddening, metal- 
licity and distance associated with cluster studies. This has 
allowed a determination of the cluster age with an uncer- 
tainty of only 0.5Gyr (errors due to mass, radius and metal- 
licity) , by comparing the location of the primary component 
to isochrones in the mass-radius plane. However, the abso- 
lute age(s) for the cluster derived from the Victoria-Regina, 
BaSTI and Y 2 isochrones differ by up to 1.3Gyr, and yet 
the two extreme isochrones overlap almost perfectly in a 
M — R diagram (Fig. [T2"|) . Until these discrepancies are re- 
solved, this prevents us from drawing a firm conclusion on 
the true value for the cluster age. 

The error in the age due to the uncertainty in the mass 
and radius is only 0.3Gyr for a single isochrone set. This 
is ~ 3 times better than what is typically obtained using 
isochrone-fitting using CMDs. 

By combining the stellar radii with a temperatures es- 
tablished from the photometry, we determined the cluster 
distance to be 4 kpc, with good agreement between the 
value from each binary component. This was used to com- 
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pare the cluster CMD to isochrones, for an adopted red- 
dening of E(B — V) = 0.15. We found that for the derived 
distance and age, the three isochrone sets match the ob- 
served CMD reasonably well - however the obtainable age 
precision is not nearly as good as when using the binary 
components in a mass-radius diagram. 

Our main conclusion is that with the needed precision 
of reddening, metallicity and distance estimates, accurate 
ages of clusters, based on "isochrone-fitting" , will be very 
difficult to obtain until GAIA can provide high accuracy 
parallaxes. 

In order to improve on cluster age-determinations, many 
more eclipsing systems in clusters should be investigated 
(Kaluzny et al. I2006|) . It would be particularly interesting 
to study clusters with more than one detached system (see 
eg. Southworth et al. 200||), as this would provide more 
constraints on the shapes of the isochrones in the M — R 
diagram at fixed age. A set of well studied clusters with 
good CMDs and accurate masses and radii for several de- 
tached eclipsing systems would provide an excellent cali- 
bration set against which stellar models can be tested to 
high accuracy. By extending such studies to globular clus- 
ters detached eclipsing binaries could ultimately provide 
the strongest constraints on their ages. We have shown in 
this study that using 8m-class telescopes it is feasible to 
work at the required level of accuracy at V = 20 thus 
making the turnoff stars in the nearest globular clusters, 
and most open clusters, available. 

Finally, it is worthwhile to note, that with further ob- 
servations of V20, it will be possible to reduce the size of 
the crrorbars for both the mass and radius determination 
and thus provide even tighter constraints (for a given set of 
models) on the age of NGC6791. With currently existing 
instrumentation it should be possible to reduce the errors 
in mass and radius to below 0.5% for both the primary and 
secondary components. 
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